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Abstract 

o ■ 

£NJ ■ The genuine twist-3 quark-gluon (qGq) contributions to the Generalized Parton Dis- 

C\ ■ tributions (GPDs) are estimated in the model of the instanton vacuum. These twist-3 

effects are found to be parametrically suppressed relative to the "kinematical" twist-3 
ones due to the small packing fraction of the instanton vacuum. We derived exact sum 
rules for the twist-3 GPDs. 



> : 1 Twist-3 DVCS amplitude 

o 

1 Deeply Virtual Compton Scattering (DVCS) |UI2] is a subject of rather intensive investigations 
during the last few years (for a detailed recent review of DVCS see [3 ). It is a two-photon 
(Nj . process: 

O; 7 *( ?)+i \r(p)^ T '( g ') + A r'(P'), (1) 

r-j ■ 

in which the initial virtual photon 7*, being produced by a lepton, interacts with a parton of 
the initial nucleon N, the final real photon 7' is radiated by the parton, and the nucleon N' in 
the final state is formed by the active parton and the rest of the initial nucleon N. The reaction 
is considered in the Bjorken limit, when the virtuality of the initial photon Q 2 = —q 2 P 2 is 
large, Bjorken variable XBj = Q 2 /2(P ■ q) is fixed and the momentum transfer squared is small: 
(P' - P) 2 < Q 2 . The first experimental results on DVCS can be found in Refs. 000111. 

We will work in the infinite momentum frame defined by the light-cone vectors n M = 
1/(P +V / 2)(1, 0,0,-1) and p» = P+/y/2(l,0,0,l), with P> = (P_+ P'Y/2 being the aver- 
age hadron momentum (note, that n 2 = p 2 = and n ■ p — n ■ P — 1). Let us define the 
momentum transfer as A M = (P' — P)^ and its square as t = (P' — P) 2 . In this paper we shall 
use a kinematical variable £, analogous to the Bjorken one, defined as £ = — (n ■ A)/2. A more 
detailed description of the DVCS kinematics can be found, e.g., in 0. 

To leading order in the electromagnetic coupling constant, the amplitude of DVCS is pro- 
portional to the hadronic tensor H^, 

Hfiv = -iJ d'xe-^ (P'\TJ^x)J u (0)\P) . (2) 

To order (1/Q)°, it can be expressed through the non-diagonal hadron matrix elements T a and 
T a of gauge-invariant light-cone bilocal operators: 
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(let us note, that here we do not write explicitly the phase factors which make the bilocal 
operators to be those gauge- invariant). The coefficient functions C±(x, £) are 



± 



x — £, + ie x + ^ — is 



and the transverse tensors —g~L, and e^, are 



' hi/ 



~ 9 ay = ~9iw + n^p v + n v p^ , 
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In general, the functions T a and T a contain both the purely quark and mixed quark-gluon 
contributions. However, in the considered case, the quark-gluon contribution vanishes due to 
contraction with the light-cone vector n a . Therefore the scalar products n a T a and n a T a can 
be parameterized in terms of four functions H q (x,£,,t), E q (x,£,,t) and H q (x,£,,t), E q (x,£,,t) in 
the following way 



n a T a 
n a T n 



H q (x, f , t) N(P', S')nN(P, S) + EJx, £, t) N(P', S 



2M 



H q (x, e, t) N(P>, S')nj 5 N(P, S) + EJx, £, t) N(P', S 
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N(P,S), (8) 
N(P,S). (9) 



Here N(P, S) is the Dirac bispinor of a hadron (a proton in our case), M is the hadron mass 
and we use the notation x = x^j^ throughout the paper. The functions H, E and H, E are the 
Generalized Parton Distributions (GPD) (for a review see Refs. jKJUIJE]]). In the forward limit 
P' = P the distributions H q and H q coincide with the usual parton densities q(x) and Aq(x). 
The distributions E q and E q are completely new ones. Thus the investigation of DVCS allows 
us to get an additional information about the structure of hadrons. In particular, the GPDs 
contain an information about the hadron spin distribution among the partons of the hadron 

0HU. 

The amplitude of DVCS has corrections of two types: firstly, loop corrections in powers of 
the electromagnetic and strong coupling constants and, secondly, higher twist corrections in 
powers of 1/Q. We will not consider the loop corrections in this paper. Making allowance for 
the first correction in powers of 1/Q, the hadronic tensor takes the following form [T21 ITS] 
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(for the pion target see JEHU)]): 
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Note that satisfies the transversity conditions 

(fH^ = 0, H^q') 



0. 



(10) 



(11) 



We see that the correction contains a term proportional to the transverse part of the vector 
functions T a and T a . This means that, in comparison to the leading order, to order 1/Q we 
deal with corrections of two types: those "kinematical" associated with purely quark operators 
and "genuine twist-3" ones originated from the quark-gluon operators, which have the form 
ipGip, where ip is the operator of the quark field and G denotes the operator of the gluon field 
strength tensor. The "kinematical" contribution can be expressed through the twist-2 quark 
GPDs H, E, H, E US HUEI E3 HZJ- 

In the case of polarized DIS the quark-gluon correction is supposed to be small relative 
to the "kinematical" one. This assumption is known as Wadzura-Wilczek approximation |18j . 
The experimental data ^OIEOIE] point at its reliability. It is assumed the same approximation 
to be valid in the case of DVCS. The aim of this paper is to test this hypothesis in the model 
of the instanton vacuum j22j . 

A method based on the model of the instanton vacuum, which allows to calculate hadronic 
matrix elements of the quark-gluon operators, was suggested [21]. In the case of DIS it was 
shown [23J|2H] that the operators of twist-3 are of order (p/R) 4 ln(p/R), where p is an average 
size of instantons and R is an average distance between instantons in the instanton medium. 
Thus, the matrix elements are parametrically small due to the packing fraction p/R~ 1/3. In 
this paper we apply this analysis to the nondiagonal hadronic matrix elements of the quark- 
gluon operators in the framework of DVCS. We show that in the model of the instanton vacuum 
they are parametrically suppressed by the packing fraction of the instanton vacuum. In the 
forward limit we reproduce the results for the diagonal case, that is for polarized DIS. 



2 "Genuine twist-3" effects 

Let us consider the following vector and axial vector bilocal operators 

V a (x,-x) = $(x)j a [x, -x]if)(-x) , 
A a (x,-x) = ip(x)j a ^ 5 [x, -x)ip(-x) . 



(12) 
(13) 
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Here 



[x, y] = Pexp <^ig J du (x — y) v A u {ux + (1 — u)y) j (14) 

is the phase factor ordered along a straight line connecting the points x and y. In what follows 
we will not write it explicitly assuming its presence in all non-local gauge invariant operators. 
In oder to separate the "genuine twist-3" contributions to the operators (fT2|). (fT3j) from the 
"kinematical" ones we shall use a technique elaborated in [T3J EH EE3 UH HZ] • According to the 
result obtained by Balitsky and Braun fIG\ l2"Tj in the case of massless fermion operators there 
are the following operator identities (£0123 = 

d f 1 - 

ip(x)-y a ip(—x) = -^-^ J duip{ux)xip{—ux) 



ie r 



sxa / duux l3 D x [i/j(ux)'y a/ y 5 ip(—ux)~\ + / du(5 a (ux), (15) 



d 

i'ixha^i-x) = 7^ / duip(ux)x^f 5 l/)(-Ux) 



afi\a \ duux /3 Q x [ip(ux)'y a 'il ! {— ux )] + / du <3 5a (ux) , (16) 



where 



•1 

<&a(x) = 9 I dv^{x)[e a ^ UT G^ p {vx)xpX u Yl5 — i v G a p{vx)x p x] , il)(—x) ) (17) 



®5q(^) = 9 j dvip(x) [s af , UT G tip (vx)x p x v YlB - ivG ap {vx)x p x\^(-x) . (18) 

In these equations the derivative D x is defined for any combination T of Dirac 7-matrices as 

d 

D x {ift(x)r[x, — x\ip(— x)\ = lim — — {ijj(x + y)T[x + y, —x + y\ip{—x + y)\ . (19) 

y^o oy\ 

Later we will use the fact that the non-diagonal matrix element of (|T9"|) is proportional to the 
momentum transfer: 

(P', S' |2) A {$(x)T[x, -x]ip(-x)}\ P, S) = i(P'-P) x (P', S' \$(x)T[x, -x]ip(-x)\ P, S) . (20) 

If we took the fermion mass into account, there would appear an additional mass term in the 
expression (fTT)|) for the axial vector bilocal operator 

+ 2im / duu ijj(ux)ia a /3X l3 'y 5 'ijj(—ux). (21) 



We would now like to express the operators (fT2"j) . (fT3|) through the symmetrical ones: 

V sym (x, -x) = $(x)x[x, -x]ip(-x) , (22) 
A sym (x,-x) = ${x)x^ b [x, -x)i/;(-x) . (23) 
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These operators are "symmetrical" in a sense that on light cone (x 2 = 0) they expand in series of 
symmetrical traceless local operators. In order to do so we should solve the system of ordinary 
differential equations (fT5|) . (|16|) with respect to V a (x, —x) and A a (x, —x) (see Appendix of ref. 
[T7] for details). The vector solution can be presented in the following form: 

V a (x, -x) = V™ w (x, -x) + V* w3 (x, -x) , (24) 

where ^ 

V™ w (x, -x) = - J du j e (i-«)[(*-®) 2 -* 2 2> 2 ] 1/2 _|_ e -(i-M)[(*-2>) 2 -* 2 £ 2 ] 1/2 J 

ru d - 

\x a D 2 — (x ■ D)Da\ J dvv ip(vx)x[vx, —vx]ij>(—vx) + -g—^ [ip(ux)x[ux, — ux]ip(— ux)~\ 

— ^a/3Ao-^ /3 S) A o — / dv i()(vx)x r Y5[vx, — vx\ij)(— vx) > , (25) 
ox a J J 

V* w3 (x, -x) = - f 1 du ( e (i-«)[(^^) 2 -^ 2 s 2 ] 1/2 _|_ e -^-u)[{x-^f-x^/-2 



— i [x 2 D a — x a (x ■ D)] / dvv / dt ip(yx)[vx, tx]gG llu (tx)x fJ,/ y u [tx, — vx]ip(— vx) 

+ <5 a (ux) - ie aPXa x p ® x dv <&l{vx) j . (26) 

Here V^ w is the "kinematical" (or Wandzura-Wilczek) contribution to the bilocal vector oper- 
ator (fT2^1 . The WW-contribution for the spin-| target was discussed in details in Refs. [Tr^lTI)]. 
The second term V^ w3 in (|24|) is the "genuine twist-3" (quark-gluon) contribution to the vec- 
tor operator. In the case x 2 = the formulae ()24|) -()26 |) reproduce the result of Belitsky and 
Miiller [T3*] . It is of necessity to note that the "genuine twist-3" contribution to the scalar 
product x a V a (x, —x) vanishes, with the first term of the "kinematical" contribution only being 
survived. In the case when x a = Xn a /2 this means that the quark-gluon contribution to the 
function n a T a vanishes, giving rise to the decomposition (J5j). 

The corresponding expressions for the axial vector operator (fT3*j) can be obtained by the 
following substitution in Eqs. (J2IJ)-(|2S|) : 

r-0 — ► , 



3 Diagonal nucleon matrix elements 



Let us consider polarized structure functions defined as integrals over diagonal matrix elements 
of string operators (for a review see [23123 and references therein): 
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dx 



iXx 



2M J 2tt 
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+ h L (x)M(S ■ n)(p^n u - p u n^) 
+ h 3 (x)M(S ±IM n u - S ±u n^), (28) 

where the spin 4- vector is normalized as 5* 2 = — 1. Taking the matrix element of (fTBj) with the 
mass term (j21j) we arrive at formulas known as Wandzura-Wilczek relations (PHI ED])- For the 
region < x < 1 we have: 



S± a { gr(x 

1 f 1 dy f dX 
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(29) 



An analogous equation for the region —1 < x < reads: 

hAx) 



+oo 

2M ]_ x y J 2vr \ 



"2"; 



(30) 



It is easy to get the Wanzsura-Wilczek relations for the Mellin moments. For example, in the 
massless case they have the form: 



S \ rv I dx X 



9t{x) 



m + 1 

In particular, for m = we have: 
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y dx [g T (x) - gi{x)) = 0; 



for m = 1 we get: 



and for m = 2 we arrive at: 
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The second term in (J34J) is purely longitudinal in a and does not give any contribution, therefore 
the transverse part of (|34j) takes the form 



Si 



9t{x) - -9i{x) 
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a = 1,2. 



(35) 



The validity of (|32 j) - (}3"4*)) follows from the observation that &^(\n/2) ~ A 2 and, consequently, 
the first non- vanishing moment in (J31|) is the x 2 -moment. The eqs. (|32j) . (|33|) known as 
Burkhardt-Cottingham |3T] and Efremov-Leader-Teryaev (33 sum rules correspondingly. Using 
the definition of the constant d^ 2 \ which measures the effect of genuine twist-3 contributions: 



dec oc 



9t{x) - ^9i{x) 



we come to the identity 



P,S> 



-AMd {2 \ 



(36) 



(37) 



which illustrates that the constant d^> characterizes the correlations of quarks and gluons in a 
hadron. 



4 Nondiagonal nucleon matrix elements 



In this section we would like to concentrate on gluon part V* w3 (x, —x) of eqs. (j2*3j) - (j2T)j) with 
x = | An. Fourier transformation of its nondiagonal nucleon matrix element 



(p', s' ; p,s ]X ) = JQ e-*> (p, s' |yr (y» ~t 



p,s 



(31 



is the gluon part of the vector GPD, where the spin 4- vectors are normalized as S 2 = S' 2 = — 1. 
Mellin moments of T 1 ™ 3 are given by 



1 d 
dx x m ^\P\ S'; P, S; x) = (-i)" 



P',S' 



V, 



tw3 
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(39) 



A=0 



It is not difficult to predict the results of the calculation of the first moments in (J39j) . The 
(J5-term in ()26|) is proportional to x 2 ()17|) and, consequently, to A 2 . The first and the third terms 
in (|26|) are proportional to x 3 and, consequently, to A 3 . This means that x°- and x 1 -moments 
should vanish; the (5-term defines the x 2 -moment, with the other terms of ()26|) being given 
their contributions to the highest moments only. Thus one can get by direct calculation that 



J dxT^ix- P', S'; P,S) = j dxxFi w3 (x; P', S'; P, S) = 0, 



(40) 



These results mean that the first two moments of the GPDs do not get gluon contributions, 
the WW approximation turns out to be exact for them. 
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The first nontrivial Mellin moment sensitive to the genuine twist-3 contributions is the 
x 2 -moment. It has the following form: 



dxx 2 ^ a w3 (x;P',S';P,S) = \(P'.S' 

1 



The transverse part of eq. (}4~Tj) has the form 



G n anP 



pnrn 



- n a G X f3n f3 -f X 



dxx 2 ^ Q w3 (x; P', S'; P, S) = - (p\ S' g^G^nPn^ 
1 x 
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(41) 



a = 1,2. (42) 



The expressions for the axial moments can be obtained using the substitution: r^i — > T^ip. 



5 Estimate of the nucleon matrix elements 

In order to estimate the magnitude of the second moments (J3Hj) . P2*|L we should analyze the 
nondiagonal hadron matrix elements of the quark-gluon operators: 



0± = gtpG ±x n x nip, 
5± = gi)G LX n x n^il) 



(43) 
(44) 



(in this section we will use a notation: 0± = O a , a = 1,2). For this purpose we first consider 
the nondiagonal quark matrix elements of the operators (|43j) . (|44jl . A method of calculation of 
such matrix elements in the model of the instanton vacuum was suggested in Refs. 
It was developed in the framework of the Effective Chiral Theory of the nucleon [33] (see also 
recent reviews |34[ l3*o] and the references therein). The essence of the method is in the following. 

In the model of the instanton vacuum [221 the fermion partition function in the presence 
of iV_|_ instantons and N_ anti-instantons can be represented as a functional integral over the 
fermion fields with a fermion effective action S e ^\^,ip\: 



7fermions 



Z N± C / VipVipexp {iS eS [ip, , 



where in the case of one quark flavor 



d A k 

W) 4 



ip{k) k-M(ik) 7p(k) 



(45) 



(46) 



and Zn ± is the variational partition function of the gluodynamics. Eq. (J4T?|) points out that a 
quark propagating in the instanton medium acquires a dynamical momentum-dependent mass 
M(ik). It is this mass that leads to the nonvanishing value of the vacuum condensate (ipip) and, 
consequently, to the spontaneous chiral symmetry breaking in QCD [31]. The quark propagator 
in the instanton vacuum model can be written in the momentum representation as 



S(k) 



k + M(ik) 
' k 2 - M 2 {ik) + ie' 



(47) 
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The dynamical mass M(ik) has the form: 



M(ik) = M\F(ik)\ 2 . 



(4* 



Here the mass M ~ p/R 2 ~ 345 MeV, where p is the average size of the instantons in the 
instanton medium and R is the average distance between instantons. The instanton medium is 
characterized by its "packing fraction" p/R with a phenomenological value p/R ~ 1/3. Thus, 
the product Mp ~ p 2 /R 2 is parametrically small. The function F(k) in ()48|) is a form factor 
proportional to the Fourier transformation of the wave function of the fermion zero-mode: 



F(k) 
F(k) 



-tj t [I Q {t)K,{t)-h(t)Ki{t)] 
1, p^O, 



(49) 



where I n (t) are the modified Bessel functions and K n (t) are the modified Hankel functions. 

The effective action (J4T)|) can be used in the calculation of the quark matrix elements of the 
quark-gluon operators PHjl - fP^ . the operators being replaced by the effective ones: 



(p',s'\0\p,s) 



inst.vac. 



tf,s'\ u O n \p,s) 



eff • 



(50) 



Here "O" is the fermion operator, which represents the quark-gluon operator O in the effective 
fermion theory. The explicit form of "O" was obtained in (2S1 EHj • The effective action and 
effective operators being known, it is possible to calculate the quark matrix elements of interest. 

The direct calculation for the case of the operators ()43|) . ()44|) results in the following ex- 
pressions: 



(p',s'\"0 ± »\p,s) 



1 



off 



Mp 2 



d A k 



G{ik) 



u(p', s') F*{i{k — p))F{ip) h^ 5 S{k — p)T*±(k, n) 

— F*{ip')F(i(k — p')) Sj_(fc, n)S{k — p')n^ 5 u(p, s) 
+ (corrections), 



(51) 



{p',s'\«0 5± n \p,s) 



off 



d A k 



G(ik) 



u(p', s') F*(i(k — p))F{ip) nS{k — p)T,±(k, n) 

- F*(ip')F(i(k - p')) ± ± (k, n)S(k - p')n\ u{p, s) 
+ (corrections). 



(52) 



Here G(k) is the form factor proportional to the Fourier transformation of the one-instanton 
dual field strength 



G(k) 
G{k) 



32tH 



-4tH 



4 



K (t) 
0. 
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8 



t = kp, 



(53) 
(54) 
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and the following definitions have been used: 

E a (A;,7ij = —gr n <V + ~^-n x a au - -n a au , 

i 1 
<V = ^bi^lA, n x = - — = (1; 0,0,-1), 
1 P+V 2 

u(p, s) - is the quark Dirac bispinor. 

Using equations of motion, the expressions (jKTj) . (JH2J) can be reduced to the form: 



(p,s \"0± n \p, s) ea = MI(p) u(p ,s)nj±'y 5 u(p,s), (55) 
(j>',s'\ u 5± "\ P ,s) eS = -^MI(p)u(p , ,s')h 1± u(p,s), (56) 



where £ = — |(p' — p) -n, and the integral /(p) is defined in leading oder of p /R as 



J(p) = p 



2 



1 /" d 4 k G(k)F 3 (k-p) 



3 J (2tt) 4 (k-p) 2 + M 2 F 2 (k-p) 



- 1+4 



(A; • p) 
k 2 p 2 



(57) 



where the integration is performed in Euclidean space with p 2 = —M 2 . In the case of diagonal 
matrix elements we reproduce (up to a sign) the results of Ref. [24J: 

{p,s\«0 ± »\p,s) eS = 4M4 2 u } arkS± , (58) 

(p,s|"0 5± "|p, S ) eff = 0, (59) 

where 

41 k = \m (so) 

is a quark constant analogous to the hadron constant (jSEj) in Eq. (pT7jl . Thus, the nondiagonal 
quark matrix elements (|53|) . (jSE} of the twist-3 quark-gluon operators can be expressed through 
the quark constant: 

(p\s'\ u O^\p,s) eS = 2Md^l Tk u( P ',s')h 1±l5 u(p,s), (61) 
(p's'lWIP.^ = -2^Md ( X k u(p',s')n 1± u(p,s). (62) 

In the limit of a small average size of instantons p we get the following estimate of the integral 

Hp) m- 

41k ~ Hp) ~ p 2 ^ 2 lnp 2 M 2 , p « A (63) 

The quantity Mp ~ 7ip 2 /R 2 is a small parameter in the effective fermion theory and, thus, 
we notice the quark matrix elements of the quark-gluon operators to be parametrically sup- 
pressed as a forth power of the "packing fraction" p/R of the instanton vacuum. The numerical 
calculation of J(p) gives d^ ark = 0.011 at the phenomenological value Mp = 0.58. 

The non-zero value of d^ aTk appears only at the order ~ n 2 p A /R A . At the present stage 
of development of our theory we cannot make precise predictions at this order, as this would 
require to go beyond the one-instanton approximation. Therefore we consider the value of 
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^quark = 0.011 only as an order of magnitude estimate. Using this order of magnitude estimate 
we can predict the order of magnitude of for the nucleon [2%] : 

^-lO" 3 , with d® - df ~ 1/N C , d^ + d^-l/Nl (64) 

Without the suppression due to the instanton packing fraction a "natural" value of the nucleon 
d^ would be of order ~ J dx x 2 g\{x) ~ few x 10~ 2 . A recent experiment at SLAC [THj quotes 

the values df ] = (3.2 ± 1.7) ■ 10~ 3 and d ( n ] = (7.9 ±4.8) ■ 10~ 3 at Q 2 « 5 GeV 2 . These data hint 
on possible "anomalous suppression" of d^ as predicted in the model of the instanton vacuum 

Using the results ()61|) and (jB*2*|) for the quark matrix elements we can make a simple estimate 
of the corresponding nucleon matrix element at A 2 = 0. In such a kind of estimate we assume 
that the nucleon matrix element at A 2 = has the same structure as in Eqs. (|61l62j) . Such 
an assumption corresponds to the use of a simple quark model to "translate" the quark matrix 
elements to the nucleon ones. More sophisticated estimates can be done with help of the 
chiral quark-soliton model, see e.g. Refs. [23 EH]. Note, however, that for the very fact of 
the suppression of the twist-3 nucleon matrix elements the usage of the simple quark model is 
sufficient. In the simple quark model of the nucleon we arrive at: 

(P , ,S'\ u Oj_»\P,S) <£ = 2M N d® N(P , ,S r )hj±j B N(P,S), (65) 
(P',S'\«0 5± »\P,S) eS = -2£M N d® N(P',S')n lx N(P,S), (66) 

where N(P, S) is the nucleon Dirac bispinor. We come to the conclusion that the hadron matrix 
elements of the "genuine twist-3" (quark-gluon) opeators are parametrically suppressed in the 
effective chiral theory by the "packing fraction" of the instanton vacuum. 



6 GPDs decomposition 

The quark twist-2 (WW) contributions to the nucleon GPD T\_ and T± have been obtained 
in Refs. |13| ITE] . In this section we would like to compare the WW contributions with the 
genuine twist-3 ones, calculated in the previous section (to be precise, we will compare their 
x 2 -moments). Following Ref. ^2| we present the "transverse" part of the GPDs in the form 1 : 

^(s,e,A) = N(P',S'){(H + E) ^ 

+ G 1 ^ + G 2 ^ + G 3 Ain + G,ie^Aln l5 }N(P,S), (67) 
^(x,e,A) = N(P',S'){H^ l5 + E^ 5 

+ G x |g:75 + G 2 7^75 + G 3 A^ 75 + G 4 ie^n} iV(P, S). (68) 

Note that in the above equations all functions depend on x, £ and A 2 . The first three moments 
of the quark contributions have the following structure. 

1 We slightly changed the notations of Ref. an d added supplementary functions which were missed in 
Ref. [TJ], see also jUJ. 
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The x°-moments are: 



f dxF^x^A) = N(P',S'){g m (A 2 ) [^]}n(P,S), 



(69) 



dxF ±fl {x,£,A) = N(P',S'){G P (A 2 



Here the well-known sum rules HU have been used: 



2M 



75 



Ga(A 2 ) [^ 7fi ] iV(P,<?). 



(70) 



y dxH(x,£,A 2 ) 
J^dxE(x^,A 2 ) 
dxH(x,£, A 2 ) 



-l 
i 



dxE(x,£, A 2 



^i(A 2 ), 
i^ 2 (A 2 ), 
GU(A 2 ), 
Gp(A 2 ), 



(71) 



where -F\(A 2 ) and F 2 (A 2 ) are the Dirac and Pauli form factors, Gp(A 2 ) and Ga{A 2 ) are the 
pseudo-scalar and axial-vector form factors; and, by definition, the magnetic form factor is 



G M (A 2 )=F 1 (A 2 ) + F 2 (A 2 ) 
The results and (J7UJ) imply that 



(72) 



dx Gi(x, £, A) = 0, 



dx Gi(x,£,A) = 0. 



These sum rules can be considered as a non-forward generalization of Burkhardt-Cottingham 
EH sum rules 

For the x-moments of the vector GPD we get: 



dxxG\(x, £, A) 
dxxG2(x, £, A) 
dxxGs(x,!;, A) 
dxxG±(x, £, A) 



1 d_ 
2d£ 
1 



dxx E(x, £, A), 



2 
0. 

0. 



G A (A 2 )- j dxx(H + E)(x, £, A) 



(73) 
(74) 
(75) 
(76) 



Note that the x-moment of the GPD G 2 does not vanish in the forward limit and gives the 
quark orbital momentum of the nucleon |12j : 



lim 



f 1 1 

J dxxG 2 (x,£,A) = -J q + -Aq = -L q . 
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For the x-moments of the axial- vector GPD we have: 



dxxGi(x, £, A) 



dx x G 2 {x,^, A) 



dxx G 3 (x, £, A) 



dxxG±(x, £, A) 



£G E (A 2 



A 2 



dxx E(x, £, A) 



4M 2 



F 2 (A 2 



dxx H(x, £, A) 



where 



- A G E (A*) 



G E (A 2 ) = F 1 (A 2 ) + ^F 2 (A 2 ) 



(77) 
(78) 
(79) 
(80) 

(81) 



is the electric form factor. These moments have no gluon contribution due to (|40j) and they have 
already been discussed in [T2] and [T^]- Similarly to the previous case one can consider the 
relations (|73 |) -(|76 |) and (|77j) - (|5D*|) as a non-forward generalization of Efremov-Leader-Teryaev 
[H2] sum rule (|33|). 

We are passing now to a consideration of the a; 2 -moments. For these moments the "genuine 
twist-3" (quark-gluon) contribution (|32Jl does not vanish and one can decompose the GPDs G 
and G from (|B7|). (|55j) into the WW and genuine twist-3 parts: 



Gf(£,A) = G^ W (£,A) + G^(£,A), i = 1,2,3,4. 

For the a; 2 -moments of the vector GPD "kinematical" part we get: 
•i 



(2)tw3, 



(82) 



dxx 2 G? w (x,Z,A) 



dxx 2 G 2 vw (x^,A) 



dxx 2 Gf w (x,£,A) 



^f 2 (a 2 ) + e ( - 1) a) - ^ (2) (e, a) 



^M(A 2) + ^(e|-l)^(e,A) 



£,Gm{A 2 



AM 2 

+ H {l \^A)-2{H + Ef 2 \^A) 
A 2 \ 



1 



dxx 2 Gf w (x,£,A) 



4M 2 M4- i )^' A) 

-^(#+£) (2) (e,A) 

^ (1) (e,A) + ^«(e,A) 



53) 
M) 

^5) 
36) 



where we used shorthand notations for the Mellin moments of the twist- 2 vector GPDs: 



H (m) (Z,A) = J dxx m H(x,£,A), 



dx x m E(x, £, A); 
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and analogous notations for the twist-2 pseudovector GPDs H and E. 
For the x 2 -moments of the axial-vector GPD we get: 



dxx 2 G? w (x,£,A) 
dxx 2 GY w {x,£,A) 



4- 2 



*7) 



dxx 2 Gf w (x,t,A) 



-1 



e(G A (A 2 ) + (H + E)^ ( ^ A) 



A ' ^ -2 )£«(£, A) -2#^, A) 



4M 2 V«9£ 
^(G A (A a ) + (lT + £0 (1) (e,A) 
A 2 \ 



1 



4M 2 



4- 2 



^ (1) (e,A) + ^«(e,A) 



(89) 
(90) 



Let us consider the quark-gluon contributions to the moments of J r ± and J r ±. We have 
already shown that their x°- and x 1 -moments vanish (|40jh with the x 2 -moment being given by 
f|42j) in the vector case. Substituting into (jUJ) and decomposing the result according to 
f)67jl we get the following expressions for the x 2 -moments of the genuine twist-3 vector GPD at 
A 2 = 0: 



-i 



dxx 2 G\ w3 (x,£) 
dxx 2 G t 2 w3 {x,£) 
J dxx 2 Gl w \x,0 
j dxx 2 G\ w3 {x,0 



0. 



-^ (2) , 



S 2 \ 



Here d^ is given by Eq. (JHUj) . For the axial- vector GPD we get at A 2 = 0: 



dxx 2 G\ wi (x,0 
dxx 2 Gi w3 {x,£) 
dxx 2 Gf 3 (x,0 
dxx 2 Gf 3 (x,0 



0. 

2 



:i-*V 2) , 



-krf (2) , 

3 s 



1 



S 2 \ 



(91) 
(92) 
(93) 
(94) 

(95) 
(96) 
(97) 
(98) 



The expressions ()91| ) -([94| ) and f|95j l -(|98 p for the the genuine twist-3 contributions should be 
compared to the expressions (J%3*|) - (j8T)j) and (|57|l - (j9*D*jl of the WW part for the quark contributions 
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correspondingly. As we saw above, in the instanton vacuum the constant is strongly 
suppressed by the packing fraction of instantons j21]. Here we expressed the x 2 -moments of the 
genuine twist-3 GPDs in terms of the constant S 2 \ see expressions ()91 )) -([9" H) an d (|95j ) -([98 p . We 
conclude that the twist-3 vector and axial-vector GPDs are parametrically suppressed relative 
to the WW ones due to the packing fraction of the instanton vacuum: 

r (2)twZ / 2 \ 2 / 2 \ 



This indicates that in the QCD vacuum the "genuine twist-3" (quark-gluon) contributions turn 
out to be small relative to the purely quark ones. 
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